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Explicit formulas for solutions of recurrence relations for 3-loop vacuum integrals are gen- 
eralized for the n-loop case. 

Recently |lj a new approach to implement recurrence relations 0] for vacuum integrals was suggested. 
Such relations connect Feynman integrals with various degrees of their denominators. In many cases they 
provide a possibility to express an integral with given degrees of denominators as a linear combination 
of a few master integrals with some coefficient functions. The common way to evaluate these functions 
is step-by-step recurrence procedure, which demands a lot of calculations. On the other hand, the 
construction of such procedure is a serious problem even at the three-loop level. 

For vacuum three-loop integrals with one non-zero mass and various numbers of massless lines the 
corresponding algorithms were constructed in [^) . 

In the approach proposed in O these coefficient functions were calculated directly as solutions of the 
corresponding recurrence relations. In pi explicit formulas for the solutions of these relations for 3-loop 
case with arbitrary masses were obtained^ As an example, the case of integrals with four equal masses and 
two massless lines was considered, and the efficiency of this approach was demonstrated by calculations 
of the 3-loop QED vacuum polarization. 

In this work we show how to extend the general formulas for the solutions of the recurrence relations 
for the multi-loop case. Let us consider L-loop vacuum integrals with TV = L(L + l)/2 denominators 
(this number of denominators provides the possibility to represent any scalar product of loop momenta 
as linear combination of the denominators; the diagrams of practical interest which usually have less 
number of denominators, can be considered as partial cases when some degrees are equal to zero): 

in 

^5 . [m D / 2 F(3- D/2)] J J D i ■ d n 

, where D a — J2iLi X^Li A aPi ' Pj ~ ^a'm 2 and pi (i = 1, . . . , L) are loop momenta. 

Let us derive recurrence relations that result from integration by parts, by letting (d/dpi) ■ Pk act on 
I | the integrand [||: 
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D6lB(n,D) = 2^^^4' ) n,I d+ (^ 1 )^ ) (I-+ Ma )S(n,Z?), (2) 
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where lfB(m, . . . , n c , . . . , ) = B(m, . . . , n c ± 1, . . . , ), and (A come from the expansion of the pi -pj 
terms appearing in the numerator through denominators D a : 

N 

p k -pi =j2( A ~%i) ( D » + ^™ 2 )- 

a=l 

Using the relations 
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the recurrence relations (0) can be represented as 
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a=l d=l 1 = 1 



1 Supported in part by the Russian Basic Research Foundation (grant N 96-01-00654), INTAS (grant 91-1180). 
Email: baikov@theory.npi.msu.su 



2 



The differential equation corresponding to (^) by substitutions ridl d+ — > d/dxd, l d — * has a 
solution <?(a; a ) = P(x a + /U a )' D ~ L_1 ' ) / 2 , where P{x a ) is the polynomial in x a of degree L: 

N 

P(x a ) = detQ^A- 1 )^) 
0=1 

Let us consider "Laurent" coefficients of the function g(x a ): 

.f\n,D ) = J ^ ¥ j...j det^- 1 )^ (*. + m*))^- 1 ^ (4) 

where integral symbols denote -/V subsequent complex integrations with contours which will be described 
below. Acting by (^|) on (0) one gets up to the surface terms the corresponding differential operator acting 
on g(x a ), that is zero. Then, the surface terms can be removed if we choose closed or ended in infinity 
point contours (to provide the cancelation of the surface terms in the last case one can consider analytical 
continuations on D from large negative values). So (0) is a solution of relations (0), and different choices 
of contours correspond to different solutions, enumerated by the index k. 

Then, according to dimensional regularization rules, integrals (Q) are equal to zero if the number of 
positive n a is less then L. In the 3-loop case |Q the coefficient functions can be constructed as linear 
combinations of those that are equal to zero if the index from definite three-index set (" Taylor" indexes) 
is not positive. One can obtain such solutions if one chooses contours, corresponding to these indexes, as 
circles around zero. In this case these three integrations can be performed and lead to coefficient in the 
common Taylor expansion in the corresponding variables. We hope that it will be applicable as well for 
the L-loop case and the number of integrations in can be reduced from L(L + l)/2 to L. 

Finally note that solutions (Q) by construction satisfy D-shifted recurrence relations 



f k (n,D) = P(I-+ f i a )f k (n,D~2) (5) 

But if one want to construct the coefficient functions f k (n,D) for some specific set of master integrals, 
one gets linear combinations of f k (n,D) with coefficients depending on D. This dependence leads to a 
nontrivial mixing matrix in the D-shiftcd recurrence relations 



(6) 



So in this paper we construct explicit formulas (^) for the solutions of the recurrence relations for 
multi-loop vacuum integrals with arbitrary masses. In forthcoming paper we plan to generalize these 
formulas for the non-vacuum case. 
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